Aggregate impacts to country and ecoregion level

library(raster)
library(sf)
source('https://raw.githubusercontent.com/oharac/src/master/R/common.R')

source(here('common_fxns.R'))

library(animation)

reload <- FALSE

1 Summary

Using the impact category level global maps, aggregate species impacts to country level. We may also want to aggregate to marine ecoregion level.

• by count (mean, variance by year; trend over all years)
• by percent of threatened species (mean, variance by year; trend over all years)

2 Methods

2.1 Set up for calcs

2.1.1 Get IDs for regions

Pull in EEZ raster and MEOW raster, and assemble into a dataframe. Include ocean area for area-weighted averaging.

2.1.2 Define helper functions

Define a function to bring in all rasters (by year) of an impact, and place into dataframe format. From the raster filename, extract the stressor/impact and the year.

2.2 Calculate impact and trend by count and percent

2.2.1 Calculate impacts by Exclusive Economic Zone

Plot the results based on the percent of local threatened species impacted by each impact category.

2.2.2 Calculate impacts by Marine Ecoregion

Plot the results based on the percent of local threatened species impacted by each impact category.